Detecting quantum-coherent nanomechanical oscillations using the current-noise 

spectrum of a double quantum dot 
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We consider a nanomechanical resonator coupled to a double quantum dot. We demonstrate 
how the finite-frequency current-noise spectrum through the double quantum dot can be used to 
distinguish classical and quantum behavior in the nearby nano-electromechanical resonator. We also 
show how the full frequency current-noise spectrum gives important information on the combined 
double quantum dot-resonator energy spectrum. Finally, we point out regimes where the quantum 
state of the resonator becomes squeezed, and also examine the cross-correlated electron-phonon 
current-noise. 



I. INTRODUCTION 

The transduction of mechanical motion of resonators 
and cantileversii^i^i^i^i^ has become increasingly impor- 
tant with the observation of motion on the nanometer- 
scale. In particular, when the ground state energy of the 
resonant mode of the mechanical system becomes larger 
than the thermal background temperature, a quantized 
state involving millions of molecules would materialize. 
As nanoelectromechanical systems (NEMS) reach this 
regime it becomes increasingly feasible, and desirable, 
to transduce their motion by coupling it to a quantum 
degree of freedom, like spinl, charge^, or flux. However, 
the challenge of finding an appropriate measuring appa- 
ratus, one whose back-action would not destroy the frag- 
ile quantum state, has not been overcome, even if such 
devices could be cooled below the quantum limit^iiS. 

Here we propose using a quantized two-level 'meso- 
scopic transport' degree of freedom, or 'transport qubit', 
as a transducer of quanta exchange, and to identify signa- 
tures of quantum coherent coupled phenomena between 
the mechanical resonator and the transport qubit. If 
successfully observed, this would validate the existence 
of a quantized mechanical state. Here we focus on a 
capacitively-coupled double quantum dot realization for 
the transport qubit. However, our analysis applies to 
several other possible devices, such as superconducting 
single-electron transistors (SSET)^ and suspended dou- 
ble quantum dotsii, which will be described later. 



A. Probing mesoscopic transport 

It is important to note that the types of experimental 
measurement that can be made on mesoscopic transport 
systems are limited; we can measure the average rate 
of particles leaving the system (current), the correlation 
between these currents at long times (the zero-frequency 
noise), and the full Fourier transform of these correla- 
tions (full frequency noise). Over the last few years, the 
zero-frequency noise has been used with great success to 



experimentally verify coherent quantum behavior (see, 
e.g., Ref. [Ill)) and may in the future serve as an en- 
tanglement measure^'^, and perhaps even aid in realizing 
a solid state test of Bell's inequalities^'*. The full fre- 
quency noise spectrum, often more difficult to measure 
in practice, is appealing because it contains information 
about the full dynamics of the system: it reveals both co- 
herent dynamics stemming from the system Hamiltonian 
H, and incoherent dynamics from the environment. This 
makes it a powerful tool for probing solid-state quantum 
systems. 



B. Summary of our results 

Our main result here is that we show how the coupled 
quantum coherent behavior, e.g., Rabi oscillations, and 
the low-energy part of the coupled double quantum dot- 
resonator spectrum, can be observed as resonances in the 
full frequency current noise spectrum. We also analyze 
the effects of temperature and decoherence on this signal, 
and show how the transition to the classical regime can 
be monitored using our approach. 

We now proceed as follows: we first define a gen- 
eral model for a transport 'qubit' coupled to the 
quantized fundamental mechanical mode of a nano- 
electromechanical resonator. This is a well-studied model 
in various forms, and has been used to illustrate, e.g., bo- 
son steering and micromaser eflFects^ii^iii. Following this, 
we explain why current-noise measurements can contain 
signatures of quantum coherent behavior. We illustrate 
this with results from our master equation model for 
two different parameter regimes. We also identify signa- 
tures of quantum state squeezingi^ of the resonator, and 
we calculate the correlation between tunneling events in 
the transport qubit and phonons leaving the mechan- 
ical resonator. Finally, we discuss other possible ex- 
perimental realizations, such as suspended double quan- 
tum dotsiiii^ii^, spin states coupled to a magnetized res- 
onator—, and capacitively-coupled superconducting single 
electron transistors 
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trons) in the left or right dot, and the state |0) rep- 
resents the empty state (N total electrons). Note that 
the excess electron in the double dot is well separated in 
energy from the other electrons due to Coulomb block- 
ade. Alternatively, the superconducting singic-clcctron 
transistor can be defined by = |2)(2| — |0)(0|, repre- 
senting the superposition of charge states on the island. 
Even though superconducting single-electron transistors 
are three terminal devices, in certain regimes the model is 
equivalent to a double quantum dot^ (see below). The 
spin-blockade case would involve a direct coupling, via 
the magnetization of the resonator—, to the electron spin 
(jz = I T)(T I ~ I I- Hereafter we retain the double 
quantum dot basis, {\L), \R)}- 

A. Master Equation 




FIG. 1: (Color online) Schematic diagram of a double quan- 
tum dot transport qubit (pink) coupled to a mechanical res- 
onator (in blue). We assume a capacitive coupling g between 
the position of the resonator and the electron charge state of 
the electron in the dot. The double quantum dot is attached 
to two electron reservoirs in the Coulomb blockade regime, 
with tunneling rates Tl and r_R. There is a coherent tun- 
neling rate A between the two charge states, and a tunable 
energy gap e. We assume that the mechanical resonator is 
already cooled to near the quantum limit^ii^. The parameter 
g is the coupling strength between the double quantum dot 
and the mechanical resonator, whose fundamental frequency 
is LOfj. In the schematic diagram, the placement of the com- 
ponents is purely illustrative. Our model is also applicable 
to circuit QED systems, where coherent energy exchange be- 
tween a resonator and a two-level system has recently been 
observed (see, e.g., Ref [iF 



II. MODEL: TRANSPORT QUBIT COUPLED 
TO A MECHANICAL MODE 



The basic Hamiltonian for a 'transport qubit' coupled 
to a quantized mechanical resonator is as follows. 



H 



(1) 



Here uji, is the fundamental frequency of the resonator, e 
is the energy gap, or splitting, of the transport qubit 
states, and A is the coherent tunneling rate between 
the two qubit states. The bosonic operators a, de- 
stroy and create excitations in the resonator. The quasi- 
spin basis describes the two possible states in our trans- 
port qubit, and we assume that the transport process 
enters and leaves through the eigenstates of Cz. For 
example, for our case of a double quantum dot in the 
Coulomb blockade regime, az = \L){L\ — \R){R\, where 
L and R represent an excess electron (N-l-1 total elec- 



Transport, in all these cases, is in non-equilibrium (left 
to right), with a large bias applied to the device, and the 
current measurement monitors the electrons/particles 
leaving the device into the right lead/reservoir (here we 
neglect displacement-current contributions). The full 
equation of motion (master equation) for this system is 
described by a super-operator Liouvillian L that defines 
the transport of particles through the "qubit" (under the 
Born-Markov approximation), bath damping and tem- 
perature terms for the resonator, 

- L[p{t)]=-^[H,p{t)]+Lo[p{t)] (2) 
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Tl and Tr are the left/right tunnehng rates, 7b is the de- 
cay rate of vibrational quanta into the resonator thermal 
bath, and T is the temperature of the resonator thermal 
bath (hereafter we set fc = ?i = 1). p{t) is the density ma- 
trix describing the state of the resonator and the qubit. 



B. Current- noise power 

We derive the counting statistics of Eq. [5] using a 
generating- function approach (Appendix A) . Using these 
equations we can calculate the current-noise power- 



s' (w)j,_ 



dT^--[{5h{t + T),5I,m^^^ (6) 
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where SIi{t) are the current fluctuations, and t ^ oo 
implies the fluctuations are around the steady-state ex- 
pectation values. This formahsm describes: 

(i) particle transport through our effective 'qubit' (i = 
j — e, electron or particle current), 

(ii) the statistics of bunching 'vibrational phonons' lost 
to the background thermal bath of the resonator {i ^ j ~ 
6, where /t, is an effective 'bosonic' current), 

(iii) correlations between electron and phonon events 
(i ^j,i = e, j = b). 

The electron current is defined by the operator 

ie = Tr sr p{t) s\^. (7) 

Similarly, the vibrational phonon current is defined by 
the operator 

ib ^ -fb a p{t) (8) 

Even though such 'phonon statistics' are typically not 
experimentally accessible, we include them here be- 
cause of the connections of our model to circuit QED 
systemai^i^i^i^i^, where the photon statistics can be 
probed with incident microwave fields and the state of 
the pseudo-spin (qubit) by suitable detectors. Such a 
system would also be suitable for observing the cross- 
correlation measurements we present later. Also, it is 
interesting to point out that in some sense the vibra- 
tional mode of the resonator itself can be thought of as 
an "acoustic phonon" with low frequency and long wave- 
length. Thus in this manuscript, for brevity we often refer 
to the "vibrational quanta of the fundamental mode of 
the resonator" as phonons. 



III. POLES IN THE CURRENT-NOISE 
FREQUENCY SPECTRUM 

To understand why the current-noise spectrum con- 
tains direct signatures of coherent quantum behavior, we 
must consider its dependence on the superoperator L. As 
discussed by Emary et and Flindt et al21 the eigen- 
values, afc, of the superoperator L, [e.g. Eq. (2)], con- 
sist of imaginary "coherent" quantum mechanical level- 
splitting terms, originating from H, and of real "incoher- 
ent" terms, originating from background thermal baths 
and non-equilibrium tunneling events. This can be seen 
by expanding the density matrix p of the coupled system 
across the eigenstates oi H, p = '^i,jN)Oli then the 
Liouvillian L acts as 

L[p] = -iHp + ipH + Lo[p] (9) 
= -i{Xi - Xj)p + Lq[p] 

As mentioned above, since all the operators in Lq are 
real, the eigenvalues of L will consist of imaginary terms 
due to energy level splitting 

SE = (A, - A,) (10) 



and real terms from operators in Lo- 
in certain conditions^, the current-noise power can 
be expanded in terms of eigenvalues at of L and the co- 
efficients Ck of the matrix {V~^IeV)kk, where Ie is the 
current operator discussed earlier, and V are the eigen- 
vectors of L, so that 

Here, Ny is the dimension of the superoperator L. If the 
incoherent terms, those outside the commutator in the 
Liouvillian L [e.g. in Eq. (2)], are much bigger than the 
coherent energy level splitting SE (e.g., Tl^r,j SE), 
then the eigenvalues of the Liouvillian are real, and the 
quantum noise is a slowly-varying function of frequency. 
If, however, the coherent terms in Eq. (2) dominate, then 
there exist poles in the current-noise spectrum around the 
absolute value of the energy level splitting — 

uj = \SE + ir + o{r/SE)\, (12) 

giving rise to the resonant features we seek. 



IV. OBSERVATION OF QUANTUM 
COHERENCE 

To illustrate how to observe quantum signatures, we 
now investigate the above model, Eq. (2), in two regimes: 
(1) an effective Jaynes-Cummings regime (when the level 
splitting matches the resonator frequency 2A = uJb), (2) 
and an off-resonance regime (where 2 A ^ Ub). Later 
on we will look at the zero-frequency noise, and make a 
comparison to a recent experiment which measured the 
zero-frequency noise of a double quantum dot in contact 
with a many-mode phonon bath. 

Hereafter, all results are calculated using the master 
equation and noise formalism described above, with a 
bosonic cut-off appropriate for the parameter regimes be- 
ing discussed. We also discuss, where appropriate, the 
dynamics of an effective pure-state, to understand how 
the energy spectrum of H contributes to the spectral 
structure of the noise. 



A. First regime: effective Jaynes-Cummings 
Hamiltonian 

An effective Jaynes-Cummings Hamiltonian can be re- 
alized if we set e = and 2A = uib- Then, 

H — Aaj; + gazia + a^) + LObo'a. (13) 

Large A implies^^ that there is a strong overlap between 
the particle wave functions in the two states, which may 
introduce extra coupling terms with the resonator. How- 
ever, for simplicity, we assume they are negligible. 
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FIG. 2: (a) shows the three important energy gaps, {SEi, 
i = 1,2,3), in the low-level energy spectrum of the Jaynes- 
Cummings Hamiltonian (see text). The parameter g is the 
coupling between the electron and the resonator, (b, c) show, 
as contour plots, how these gaps can be observed in the 
current-noise frequency spectrum S{io)e.e/2eIe, for cub = 1, 
Fl = Tfl = 0.01, A = 0.5, 7i, = 0.05 and T = 0,LOb, for (b), 
(c), respectively. The energy gaps shown in (a) are clearly 
visible as three resonances in (b) and (c). The horizontal line 
corresponds to the physical process of an electron tunneling 
without exchanging quanta with the resonator. The top and 
bottom resonant lines are proportional to the coupling g, and 
thus represent the physical process of the electron coherently 
emitting a phonon into the resonator. As the temperature 
is increased, the visibility of the two 'Rabi peaks', which are 
signatures of coherent quantum behavior of the resonator, de- 
creases. 
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FIG. 3: (color online) Normalized current-noise frequency- 
spectrum, S{uj)c,c/'iele, versus u), for ujb = 1, Fl = F_r = 
0.01, A = 0.5, 7! = 0.05, and T = 0, 0.5L^b, Ub (figures (a), 
(b), and (c), respectively) and a selection of electron-resonator 
coupling strengths g; black, green, red are (? = 0, 0.2, 0.4, re- 
spectively. The resonances are marked by arrows. As the tem- 
perature is increased, the system enters the classical regime 
and the two Rabi peaks decrease and become hard to distin- 
guish. 



First, we write the diagonal energy term for the qubit 
(ctz) in the off-diagonal basis {ax) by substituting raising 
and lowering operators in that basis 



lay 



'-{a^+iay). 



(14) 
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Then performing the rotating-wave approximation in this 
basis, by dropping counter-rotating terms, we obtain 



form resonances in the noise frequency spectrum. In par- 
ticular, because = here. 



Hjc ~ g{<7^a + a'') + LUba^'a + AcTx (15) 

- ujijo) a -f Ada;. 
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This has the spectrum of an infinite number of non- 
interacting multiplets with eigenstates, 

|±)„ = -^(|n,l,)±|n + l,0,), (16) 

where |n) is the number state of the mechanical res- 
onator, and \0)x and \l)x are the eigenstates of cTx (i.e., 
the bonding and anti-bonding states within the double 
quantum dot). 

If we consider the zero-temperature limit and a strong 
damping of the bath, then only the lowest number states 
of the mode n = 0, 1 strongly contribute to the transport 
processes (this case is well into the quantum regime, and 
the ideal situation). This regime is feasible if the effective 
temperature of the resonator is below huji,. In this case, if 
Tl ~ 7& then the initial state of each 'round' of transport 
would be 

m = 0)) = |0, L) = i=(|0, 0.) + |0, 1.). (17) 

The second component, |0, 1^;) couples to the n — and 
n = 1 states of the mechanical resonator via the |±)„=o 
eigenstates of Hjc- The first component, |0,02;), acts 
as an 'interaction free' transport route because it is the 
ground state oi Hjc- The component |0, O^;) has a unique 
'ground state energy' £^0 = while the two |±)„=o 
eigenstates of Hjc have energies 



where 



= cjf, - 2A. 



(18) 



(19) 



Our numerical simulations in Fig. [2] and Fig. [3] show 
clearly how the energy level splittings SEi, [i — 1,2,3) 



5Ei/2 ^ E±- Eo = 2A±g 



(20) 



are the upper and lower resonance 'branches' in Figs. 
2 and 3, caused by the coherent coupling between the 
double quantum dot and the mechanical resonator, and 
6E3 = 2A is the central resonance because of coherent 
internal oscillations within the dot alone. This occurs be- 
cause of the |0, Ox) ground state ofHjc, described above, 
which only evolves in time with a phase factor Eq = — A. 

As we increase the temperature of the mechanical res- 
onator thermal bath, the upper and lower resonance 
branches gradually disappear, and the central resonance, 
determined hy lo = 2A, dominates. Increasing the tem- 
perature of the 'bath' means that the mechanical res- 
onator would be in a thermal mixture of number states; 
thus for the electron, more transport channels become 
available. This is more clearly apparent in the mag- 
nitude of the noise shown in Fig. [31 illustrating that 
by monitoring the peaks in the current-noise transport 
5e,e(w) one can, in principle, distinguish classical and 
quantum behavior. However, the observation of near 
zero-temperature oscillations is not always proof of quan- 
tum behaviou r^'^^i^^'^°i'^^'^^ as they can also be described 
by a classical model of coupled linear oscillators. For ex- 
ample, in our current-noise formulation "false signatures" 
from interactions of the qubit with nearby classical oscil- 
lators may appear in the spectrum and be mistaken for 
quantum Rabi behavior. We discuss this further in the 
next section. 

One can understand the transition to the high- 
temperature case by assuming the initial state to be 



mt = 0)) 




(21) 



which connects each multiplet in the spectrum of the 
Jaynes-Cummings Hamiltonian with its two nearest en- 
ergy levels. The subspace of the Hamiltonian connecting 
\n— 1, la;),|n. Ox) \n, Ix) and jn-l-l, O^,) is (where the basis 
here is for ax diagonal). 



if, 



n— l,n,n+l 



^ (n - l)ujb + A gV" \ 

g^/n nLUb — A 

uLJb + A gy/n + 1 

\ gy/n + 1 {n + l)uj- A J 



r 



(22) 



Then, we easily see that the probability that the left dot 
is occupied (corresponding to the probability of the su- 
perposition of bonding and antibonding states --i=(|02;) + 



\lx))), is given by, 
PL(i)-^<^C„cos 



n=0 



-gWn + l - Vn)i - 2tA" 



(23) 
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Considering both an equal superposition, C„ = 1/VN 
(but with cut-off of the sum in Pl (t) at a given N) , and 
a coherent state distribution, C„(z) = z"e^^/n!, we ob- 
serve that the oscillations in the the probability PL{t) 
collapse over time, until only small oscillations with pe- 
riod A around Pl = 0.5 remain. This is because the 
non-commensurate Rabi frequencies in Eq. |23| interfere 
destructively. For a small number of number states iV, or 
a small coherent state distribution z, there is some revival 
in Phit), but as N increases the number of revivals fall. 
This is also true if the initial state is a separable density 
matrix with the resonator state in a thermal Boltzman 
distribution, as is the case for high-temperatures. 



B. Second Regime: Off-resonant interaction 

In the previous section we showed that on-resonance, 
2 A = Wb, the lowest part of the energy spectrum of the 
coupled system was visible in the current-noise. We can 
now verify that these resonances really stem from the low- 
energy spectrum of the Jaynes-Cummings Hamiltonian, 
and indicate coherent quantum dynamics, by inspecting 
the off-resonant legiine, 2 A ^ w^, where the energy levels 
have a hyperbolic behavior. In terms of the double-dot 
realization, we point out that assuming a small A implies 
a tight confinement of the electron within each dot. 

Observing Fig. 4(a), we can see upper and lower res- 
onance branches, but in this case (2 A ^ aj^) they have 
the typical hyperbolic tails of an avoided level crossing. 
In addition. Fig. IHb) shows that, as the coupling to the 
mechanical resonator g is increased, the gap in the level 
crossing increases. Once more we are successfully observ- 
ing the low-energy spectrum of the Hamiltonian in the 
power spectrum of the current-noise. For example, the 
upper and lower branches are simply given by the lowest 
eigenvalues of the Jaynes-Cummings Hamiltonian, 



5Ei,2 = |tJf,/2± + V/2 + A|. 

Furthermore, we note that there is an energy gap which 
halts the electron current in the limit when e = and 
when the coherent tunneling within the dots is small rel- 
ative to the coupling to the mode, A <C g. This occurs 
because the tunneling of an electron requires an energy 
loss proportional to the displacement of the mode, and 
because the rotating wave approximation is no longer 
valid. For transport to occur, the electron must tunnel 
from the left to the right state, which is now shifted in 
(relative) energy by 2g{a + a^). This becomes more and 
more difficult as the coupling g is increased, resulting in 
a "current blockade" effect. 

Finally, as discussed in the previous section, we 
point out that oscillations alone may not provide suffi- 
cient proof of quantum behavior. Recent circuit-QED 
experiments^!^ have focused on the idea of observing 
the square-root dependence of the energy of the Jaynes- 
Cummings system on the photon occupation number n, 
which is sufhciently distinct from the behavior seen in 
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FIG. 4: Current-noise frequency-spectrum, S{ui)e,e/2ele, 
versus both the tunneling rate A and frequency ui, for (a) 
g = 0.1 and (b) g = 0.4. In both cases uib = 1, Tl = = 
0.01, 76 = 0.05, e = 0, and T = 0. This was obtained nu- 
merically by solving the master equation in Eq. (3). As A 
approaches LUb/2 = 0.5 we see the three resonance points pre- 
viously shown in Fig. 3. Increasing g increases the gap be- 
tween the resonant peaks. The hyperbolic behavior comes 
from the well know Jaynes-Cummings eigenvalue spectrum, 
as recently observed experimentally in a similar system (e.g., 
Refs PSiH). 
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classical models. However, the preparation of arbitrary 
Fock states in a nano-mechanical resonator is not readily 
realizable at this point in time. 



C. Zero frequency noise: comparing the single and 
many-mode cases 

In the previous sections we showed how the low- 
energy levels of a Jaynes-Cummings Hamiltonian can be 
seen in the /uH-frequency current-noise spectrum. How- 
ever, most recent experiments have focused on the zero- 
frequency noise. For example, Kiefilich et all^ showed, by 
comparing experiment and theory, that coherent oscilla- 
tions in a double quantum dot produced sMper-Poissonian 
[S'(0)e,e/2e/e > 1] signatures in the zero-frequency noise, 
while incoherent transitions (sequential tunneling in- 
duced by increasing the temperature of the phonon bath) 
produce su6-Poissonian noise, [S'(0)e,e/2e/e < 1]. 

Mimicking their parameter regime, i.e. considering 
their device as coupled to a resonator (or phonon cav- 
ity), now we also look at the zero-frequency noise (as a 
function of double quantum dot level detuning e) . We ob- 
serve similar signatures to theirs in the noise spectrum, 
but with a more complicated structure. We also observe, 
in Figinja), that increasing the temperature of the sin- 
gle mode resonator decreases the zero frequency current- 
noise, eventually resulting in sub-Poissonian behavior. 

Similarly, increasing the bare coupling strength g to 
the single mode resonator has a drastic effect. As 
Fig. ISUb) shows, the noise profile quickly becomes sub- 
Poissonian, developing a new peak structure around e = 
1. Interestingly, the e = point, where we earlier probed 
for coherent signatures, remains around 5(0)e,e/2e/e = 
1, indicating that coherent transport is still occurring. 

Figure \M^c) illustrates the (non-normalized) cross- 
correlated noise, i.e. the correlation between electron 
tunneling events and phonons leaving the mechanical res- 
onator into a heat bath (with rate 7). As expected, there 
is no correlation between tunneling events when the sys- 
tems are uncoupled. Furthermore, the correlated noise 
is large when e = fcwf,, where k is an integer. While the 
correlated noise grows for larger 'fc', the current itself be- 
comes smaller^^. This is simply because as e increases, 
the current can only flow through phonon assisted tun- 
neling, which happens at integer numbers of the phonon 
frequency. 



V. SQUEEZING THE QUANTUM STATE OF 
THE RESONATOR 

We have shown that the electron current-noise, 
S'e,e(w), serves as a detector of coherent interactions be- 
tween the double quantum dot and the single mode of 
the mechanical resonator. Already this is a significant 
step, as 5*6,6 (w), serves as a tool for experimental observa- 
tion. However, we can proceed a step further, and briefly 
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FIG. 5: (Color Online) (a) Current-noise frequency-spectrum 
S'(0)e,e/2eJe versus e, for lj = 1, Fl = 0.1, Vr = 0.001, A = 
0.1, 76 = 0.01, g — 0.0008, increasing T from to 2 in steps of 
0.5. Increasing the temperature reduces the super-Poissonian 
character of the current-noise, (b) S{0)e,e/'2ele for uj — 1, 
Fl = 0.1, Fr = 0.001, A = 0.1, 7i, = 0.01, T = 0, for g from 
to 0.4 in steps of 0.1 (besides g = 0.3, which has been omitted 
for figure clarity) . For strong coupling we see that the current- 
noise becomes almost entirely sub-Poissonian (as indicated 
by the shaded pink region), (c) The cross-correlation S{0)e,b 
against e for the same parameters as (b). Recall that S{0)e,b 
is defined as the correlation between the electron and phonon 
currents. For g = and for negative e there is zero correlation 
between phonon and electron tunneling, as expected. 
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consider the statistics of the phonons in the mechanical 
resonator. In such a mechanical system, these quantities 
are difficult, if not impossible, to access. However, it is 
informative to understand how the phonon statistics of 
the resonator change as we increase the temperature, and 
leave the quantum regime. 



A. Squeezing signatures 

In the proposal by Rodrigues et al^ they show that 
the resonator can exhibit properties akin to a micro- 
maser, due to the nonlinear coupling to an SSET. In 
their case, the qubit is represented by a superposition 
of island charge states cr^ = |2)(2| — |0)(0|. However, 
they focused on the regime where Wb/F — 1, observing 
that this is where the interaction between the resonator 
and SSET is maximized. In the results we have shown in 
the previous sections, we assume that the quantum dots 
and leads are weakly coupled, Wb/F ^ 1. Furthermore, 
we assume that the resonator is strongly damped (e.g., 
via cooling by another SSET, or by the double quantum 
dot itselJiSii^^) , so that only the few lowest bosonic levels 
are excited. 

However, even for our 'slow' regime, we see sub- 
Poissonian signatures in the boson emission noise spec- 
trum emitted into its nearby heat bath Sb,b{(i)/'2Ib, as 
well as in the Fano factor Fq of the number state occu- 
pation n of the resonatori^^i^i^i^l 



in') 



(24) 



as shown in Fig. 6(a) and (b). Both 'measures' identify 
similar regions of squeezing, though there is a conceptual 
difference between the squeezing of the phonons emit- 
ted (dynamically) into the heat bath, and a direct mea- 
surement of the static steady-state phonon occupation 
number. Furthermore, we see that as the temperature is 
increased, both quantities increase non-linearly in mag- 
nitude. 

In addition, we consider the correlated electron-phonon 
noise. We naively expect that stronger correlations will 
occur in the quantum regime. Figure 6(c) verifies this, 
and shows that a maximum in the correlated noise occurs 
around g = 0.1, and an increase in temperature reduces 
the overall magnitude. This is an indication, continu- 
ing from previous suggestive resultsi^, that the quantum 
noise correlation between two open systems could serve 
as a measure of entanglement, though a direct correspon- 
dence has yet to be identified. 



B. Quadrature versus number-state squeezing 

The squeezing in Fig. 6 is number state squeezing, 
and a sub-Poissonian variance in n (Fq) implies anti- 
bunching of the phonon statistics^. This is only one 



of several types of squeezing. For example, in quan- 
tum optics, generalized quadrature squeezing is often in- 
vestigated. Typically the axis of squeezing might not 
been known, so a homodyne measurement of the oc- 
cupation statistics must be performed. A homodyne 
measurement^^,, using a local oscillator to introduce a rel- 
ative phase, reveals the variance of any desired quadra- 
ture. Thus, in principle, it is possible to measure the 
normal ordered squeezing via 

(: (AQf :> ^ (: :) - (Q)^ (25) 

where Q is the quadrature defined by a desired angle 4>, 
so that 

Q^ae-''f' + a^e''l'. (26) 

Squeezing of the quadrature is implied when 

(: (AQ)2 :) < (27) 

for some given (/), because of the normal ordering. Again, 
in a nanomechanical system such a measurement is not 
feasible, but has been proposed in transmission line 
resonators^. Is is trivial to see 

(: (AQ)2 :) = (a^')e^'^ + (a^)e-^"^ + 2{a^a) 

- {a^fe^'^ + {afe-^'^ + 2(a) (a^). (28) 

However, for our model and parameter space, we were 
not able to observe any instance of quadrature squeez- 
ing. In the previous sections we discussed how strong 
contributions to the steady-state solution of the mas- 
ter equation arise from the low-level Jaynes-Cummings 
eigenstates. Our results illustrate that, in our system, 
these states only produce number state squeezing in the 
resonator mode, but not quadrature squeezing?^,. 



VI. REALIZATIONS 

As mentioned before, our model can correspond to 
charge states in a double quantum dot in a capacitively- 
coupled or suspended geometry. For the suspended 
geometrjiiiii^, it has been shown that there is a direct 
coupling between the electron wave function and a sin- 
gle phonon mode because of van-Hove singularities in the 
density of states. However such experiments have not yet 
been performed in the energy regime of the fundamen- 
tal vibrational mode of a mechanical resonator. Also, 
our model is related to that of a superconducting single 
electron transistor (SSET) capacitively coupled to the 
resonator^'^'^°. Typically there are some differences in 
the transport properties as an SSET is a three-terminal 
device, and the SSET drives the resonator into complex 
types of limit-cycle behavioriML. 



A. Energy scales 

To check the feasibility of our results we need to verify 
the appropriate energy scales in real systems. We assume 
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that our state-of-the-art resonator has a fundamental fre- 
quency of = 1 GHz. The corresponding 'resonant' 
bias, e — huj, is approximately 4 ^eV. We assume we are 
near the quantum limit, i.e., kT « hujb, T « 50 mK. 
Normal capacitivc coupling strengths for an SSET are 
100 MHz, corresponding to g = uj/10. In Figure 2 we saw 
signatures of quantum coherent oscillations for this range 
of coupling strengths. The same range {g ~ w/10) is fea- 
sible for the coupling between a double quantum dot and 
the resonator (with capacitive couplingi^). The achiev- 
able coupling strengths for suspended geometries are not 
precisely known now, but because of van-Hove singulari- 
ties in the density of states one can expect large effective 
coupling strength o^^'^^ . Finally, the inter-state tunnel- 
ing, denoted by A in our discussion, is typically tunable 
for double quantum dots. Thus a range of A ~ (1~10) 
fieV is feasible. 

B. Magnetized resonator interacting with electron 

spins 

A recent proposal^ focused on a magnetized resonator 
which interacts with one of two electron spins in a spin- 
blockaded double quantum dot system. In this case, the 
current is used to measure the spin state because, if the 
two spins are parallel, current cannot flow. An oscillating 
magnetic field, from the magnetized resonator, couples to 
one of the spin states, and thus this spin plays the role of 
a 'transport qubit' in our earlier language. The question 
of cooling such a magnetized resonator and then coupling 
it to a nearby electron spin via its quantized motion, and 
henceforth the quantized magnetic field motion, has not 
been addressed. In that case, the Hamiltonian of the spin 
and the resonator is, 

Hq = - ^ + hLJba^a + C\ {a + a'')a^ (29) 

2 V 2mcsUJb 

where C = 0.16 mT/nm. This (Eq. ^ differs from the 
Hamiltonians in Eqs. [T3l and [T5l in that Eq. [29] is diago- 
nal in the qubit energy basis. The ground state motion 
of a 1 GHz resonator is 2 x 10""'^'' m, which, using the 
parameters from Ref. 0, would generate a field of just 
3.2 X 10~^ mT, a Rabi frequency of about 100 Hz, which 
is negligible in comparison to nuclear hypcrfine and spin- 
orbit effects. Optimizing device design can increase this 
Rabi frequency considerably. For example, a larger mag- 
netization could be achieved by using a Dysprosium (Dy) 
micromagnet instead of Cobalt (Co) (giving a factor of 
about two). Similarly, a larger micromagnet thickness 
could also contribute a factor of about two to the field 
felt by the electron spin. Decreasing the distance between 
the dot and resonator could contribute up to a factor of 
ten, and using a slower frequency resonator, for a larger 
ground state displacement, could add a factor of about 
five. Taking these factors into consideration gives a Rabi 
frequency in the range 10-100 kHz. This Rabi frequency 
is still, in comparison to the charge-based quantum dot 



and SSET systems, a weak coupling, and is vulnerable 
to dephasing from nuclear hyperfine fields. However, the 
future evolution of this technology may make such an ap- 
proach feasible and desirable, especially considering the 
possible benefits of combining spintronics and nanome- 
chanics. 



VII. CONCLUSIONS 

We have illustrated how quantum coherent behavior 
and the energy spectrum of a nanomechanical resonator 
can be identified using full-frequency current-noise mea- 
surements through a nearby transport qubit. In the zero- 
frequency limit, we showed that a single-mode 'environ- 
ment', as represented by a nanomechanical resonator, 
produces unique signatures that differ from those ob- 
served in multi-mode environments. Furthermore, we 
identified regimes where phonon squeezing and cross- 
correlated noise, indications of complex quantum phe- 
nomena, could occur. All of these features could be 
realized with a double quantum dot or superconduct- 
ing single-electron transistor operating as the transport 
qubit. In a broader context, we expect that noise mea- 
surements could also be useful in two-resonator circuit 
QED systema^^i^i^, which may offer an interesting area 
for future investigation. 
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APPENDIX A: NOISE FORMALISM 

To calculate the quantum noise^^ of a system with 
Hamiltonian H, and corresponding transport environ- 
ment described by a Liouvillian L, we employ a gener- 
ating function approach. The Master equation for the 
matrix elements of the generating function g is 

^g(5l, ■■•,5m,0 = M{si, ...,S,n)g{si, ...,Sjn,t), (Al) 
Ot 

which can be formally solved by diagonalizing 

M(si,...,s™) = (A2) 

V{si, Sm)D{si, Sm)V~'^{si, Sm)- 

Here M is the Liouvillian L recast as a function of the 
counting variables (si,...,Sm)- Each Sj is a continous 
variable which tracks the passage of the current through 
system i. This gives a general formalism for calculating 
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the generating function of m coherent and interacting 
transport systems, each with a single 'one-way' current 
flow. 

The next step is to use the MacDonald formula^ for 
the symmetrized noise power correlator between systems 
i and j 



dre^-^[{SI,{t + T),SIjmt^^ (A3) 



which can be written as (s — {si, S2, s„i}) 



= {ds.,s,+k]ds,) I dTsin(wT) (A4) 
X ATrG(sr)U^i, 



where an omitted term 2t(/)^ in the integral does not 
contribute in the final result obtained upon performing 
the Laplace transformation. Noting that G{s,t = 0) = 
p{0), where the initial condition p{0) is the steady state 
density matrix and using 

^G(s,r) = A/(s)G(s,t) = M{s)e^^'^'^G{s,T^ 0), 
and the spectral decomposition of M{s), one obtains 



uj''D{s) 



D{s)- 



(A5) 



where the notation \\(xi^j^,Xir^jr^, = '^i^o^u takes 
into account the trace in Eq. (|A4[) . Note that the first 
derivative in the single system correlator ds yields 2e(/), 
and therefore 9^ provides the deviation from the shot 
noise. Using the Ramo-Shockley theorem^i, the displace- 
ment current contribution can either be omitted (by as- 
suming that the capacitances of the devices are extremely 
asymmetric, so that clc_r <C 1), or calculated using a 
multi-variable approach, because the total current fiuc- 
tuations can be written as 

SI{t + T)6I{t) = a^SlL{t + T)SlLit + T) 

+ P^SlRit + T)5lRit) 

+ al3iSlL{t + T)SlR{t) 

+ SlR{t + r)SlL{t)). (A6) 

The left and right correlations are trivially calculated us- 
ing separate counting variables for each lead. 

Equation (|A5|1 allows one to calculate the noise spec- 
trum for transport through an arbitrarily complex quan- 
tum system. This can be evaluated either using finite dif- 
ference derivatives around s = 1, or following the meth- 
ods employed by Flindt et alJ&i^. In the latter case 



we can use their approach to show that, in general, the 
cross-correlator can be written as 



dt{n,it)n,it)) = Tt[L, J2 "j*"^"' 



),(«2), 



+ Tt[Lj ^ (A7) 

ni,7i2,... 

Furthermore the terms 

^ n,p("^)'("^)'-=9,^G(s,r)U.i (A8) 

ni ,712,... 

can be evaluated by Laplace transforming the equation 
of motion 

drG{s, t) = {Lo + J2 s^L^)G{s, t) (A9) 

i 

and taking derivatives in the counting variables Si, giving 



ds,G{s, -iijj)\s=i = F{~-iuj)UF{-iuj)p{Q) (AlO) 



where 



F{-iuj) = {-iuj- L)-^ 



(All) 



and p(0) is the steady-state initial condition. As shown 
by Flindt et al^ one can evaluate this inverse by writing 



where 



F{-iuj) = -P/iLU - R{uj), 
R{cu) = Q{tuj + L)-^Q, 



P = p{0) ® 1, Q = l-P. 



(A12) 
(A13) 



(AM) 



Inserting all these expressions into the cross-correlator, 
and using Pp{0) — p{0) and Qp{0) — 0, gives the noise 
power as the trace of an inverse, 



^^"^^''^ - Re{-TT[L,R{Lu)Ljp{0)]-TT[LjR{Lo)L,p{0)]} 



2e2 " 

+ (5,,,Tr[L,p(0)]. 



(A15) 



All of the above allows us to calculate the full frequency 
spectrum for an arbitrary number of coupled systems. 
In addition, it allows us to calculate phonon current and 
statistics. We choose as the phonon current operator the 
operator which absorbs a phonon number state from the 
mode and puts it in the background bath. 
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FIG. 6: (Color online) (a) Bosonic current-noise frequency- 
spectrum, Sb.b{'^)/2Ib, versus g, for uib = 1, r_L = Fr — 0.01, 
A = 0.5, 7b — 0.05, e = 0.0, and a range of temperatures. The 
pink regions [also denoted by "Sub-Poissonian regime" in (a)] 
indicate the regimes where quantum state squeezing occurs, 
for low temperature and intermediate couplings, (b) shows 
the number-state Fano factor ((n^) — {n)^)/{n) of the bosonic 
system. The zero-temperature case closely corresponds to the 
'phonon current-noise' in (a), (c) shows the electron-phonon 
correlated noise Se,b{0) versus g. Interestingly, increasing the 
temperature decreases the correlated noise. Moreover, and 
as expected, the zero coupling point {g — 0) remains around 
S{0)e,b = for all values of temperature. Also, in a very small 
regime of weak coupling {g — > 0), S{0)e.b can be negative. 



